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We investigate general structure of Landau free energy for stabilizing a novel superconducting state with both 
even-frequency and odd-frequency components in gap function. On the basis of the Luttinger-Ward functional, 
we elucidate an emergent mixing between different "parity" in time. The simplest case of the conventional 
i-wave singlet mixed with the odd-frequency triplet state under broken time-reversal symmetry is examined 
to demonstrate the anomalous structure of the free-energy functional. The induced odd-frequency component 
alters behaviors of physical quantities from those obtained by neglecting the odd-frequency component. The 
novel mixed state may also be relevant to strong-coupling superconductivity coexisting with ferromagnetism. 

PACS numbers: 74.20.-z, 74.20.Rp, 74.25.-q 



I. INTRODUCTION 

A symmetry of order parameter is one of central issues 
in study of superconductivity. The isotropic s-wave sin- 
glet pairing in the original BCS theory has been extended 
widely to anisotropic ones including spin-triplet states in the 
strongly correlated systems 1 . Alternative direction of exten- 
sion was proposed by Berezinskii 2 , where the gap function is 
"anisotropic" in time, i.e., with an odd-frequency dependence. 
This novel state has been investigated extensively in the con- 
text of the disordered Fermi liquid 3 ' 4 , the high-r c cuprates 5-7 , 
the doped triangular antiferromagnets 8 , one-dimensional or- 
ganic systems 910 , quantum critical spin fluctuations in heavy 
fermions 11 , strong-coupling local electron-phonon systems 12 , 
and orbitally degenerate systems with full spin polarization 13 . 

Besides exploration of possible odd-frequency pairings in 
homogeneous bulk systems, it has been argued that the odd- 
frequency pair amplitude arises in the spatially nonuniform 
situation, particularly at the surface/interface of the sample, 
quite ubiquitousl y 14 ' 15 . It is more promising to induce the 
odd-frequency pair amplitude in the nonuniform system with 
lower symmetry. The case of spatially lower symmetry has 
also been investigated in noncentrosymmetric systems with- 
out inversion symmetry, where even- and odd-parity compo- 
nents are mixed, and it influences many of superconducting 
properties in a profound way 16 . Quite analogously a mixing 
of even- and odd-frequency components should occur in prin- 
ciple, when the time-reversal symmetry is broken. The present 
authors proposed the mixing of the conventional s-wave sin- 
glet with the odd-frequency triplet component under applied 
uniform magnetic fields through Zeeman splitting of the con- 
duction band 17 . Such a situation may also be relevant to su- 
perconductivity coexisting with ferromagnetis m 13 ' 18 . 

When the odd-frequency component is involved, a delicate 
treatment of the superconducting state should be required. 
Namely, in describing a superconducting state, there exist 
two independent gap functions, A(fc, ia>„) and its particle-hole 



converted counterpart A + (k,ioj n ), which are usually consid- 
ered as hermite conjugate pairs. However, it has been argued 
that the relation between them is nontrivial in the case of the 
odd-frequency pairin g 19,20 . A guiding principle of determin- 
ing the relation between them is real and minimum condi- 
tion of the free energy, and it is concluded that the correct 
choice is A + (k, ia>„) = A(fc, ia> n )*^^, although it has long 
been believed A + (k,iu) n ) = -A(k,iu> n y in the case of the 
odd-frequency pairing. The latter relation was the main source 
of fatal deficiency of the bulk odd-frequency pairing, such as 
thermodynamic instability and unphysical negative Meissner 
kernel in the superconducting phasei^i 2 !. 

Nevertheless, we shall discuss in this paper that the above 
relation, A + (k, ia>„) = A{k, ioj n )*, is no longer satisfied si- 
multaneously for both components in a coexistence case of 
the even-frequency and the odd-frequency pairings 17 . Conse- 
quently, the majority component of them takes the ordinary 
sign, A + (k,iu>„) = +A(fe, *'&>«)*, while the minority compo- 
nent must have the opposite sign, A + (fc, ia>„) - -A(fe, iu) n )*, 
in order to gain interference energy between two components. 
As a result, it exhibits different behaviors in various physi- 
cal quantities from those obtained by neglecting this minority 
component with the anomalous negative sign relation. 

The paper is organized as follows. In §2, we discuss general 
structure of the free-energy functional in the coexistence case 
on the basis of the Luttinger-Ward functional theory. The real 
and minimum condition of the free energy requires the op- 
posite sign relations for spin singlet and triplet components. 
In §3, we examine the s-wave singlet with the induced odd- 
frequency triplet under the applied magnetic fields, based on 
the ordinary electron-phonon model with a single Einstein fre- 
quency. The anomalous feature of the free-energy functional 
is demonstrated in the strong-coupling limit with single gap 
approximation. We show the comparisons of the specific heat 
and the superfluid density with or without the odd-frequency 
component by solving the full Eliashberg equations. The 
last section summarizes the paper. Throughout the paper, we 
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have neglected the orbital effect of the magnetic fields. The 
derivation of the explicit form of the free-energy functional is 
given in Appendix A. The noncentrosymmetric case with the 
Rashba spin-orbit coupling is briefly discussed in Appendix 
B. 



H. GENERAL ARGUMENT ON EVEN-FREQUENCY AND 
ODD-FREQUENCY MIXING 

A. Free energy in terms of Luttinger- Ward functional 

Let us begin with the Luttinger- Ward functional for the 
thermodynamic potential 22,23 measured from non-interacting 
one, 



n[G] =- T VTr [ln[-G- l (k)}-[l-G Q l (k)G(k)}] 

k 

<£[G] 



Oo, (1) 



where k = (fc, iu>„) is the 4-dimensional momentum with the 
fermionic Matsubara frequency, a> n = (2n + l)nT at tempera- 
ture, T. The trace is taken over the 4x4 Nambu space defined 
as = [cj(T), c\(k), c^(-k), c^(— k)] (a hat represents a 

4x4 matrix). G (k) and Q = -(772) Zk Trlnt-G- 1 ^)) are 
the non-interacting Green's function and its free energy (ther- 
modynamic potential), respectively. Here, the thermodynamic 
potential is expressed in terms of the full Green's function, 
G(k). It can be shown to be stationary at the physical G(k) 
that satisfies the Dyson equation, G~ 1 (k) — G Q 1 (k) - £(&), pro- 
vided that the self energy is given by the functional derivative 
as X,y(£) = T~ l 5$>l5Gji{k). The functional O consists of so 
called the skeleton diagram of G(k) with interaction lines. 
To be specific, we consider, 



$[G] = T — £ V(k - k')Tr[G(k)p 3 G(k')p 3 ] . 



(2) 



where p 3 is the z-component of the Pauli matrix acting on 
the particle-hole space, and V(q) is an arbitrary g-dependent 
scalar interaction (V > represents an attraction), which is a 
real and even function of q. It can easily be shown that the 
stationary condition gives the Eliashberg equation 24 with the 
self energy, 



£(k) = T^V(k-k')p 3 G(k')p 3 



(3) 



In order to perform the following Landau free-energy ar- 
gument, it is more convenient to transform Q.[G] to a func- 
tional with respect to the self energy, Q[£], It can be accom- 
plished by the Legendre transformation 25 , c£>[G] = ©[£] + 
T Z k Tr[G(k)£(_k)] as 

= ~ J] Tr In [- \G Q l {k) - t{k)}\ + ^ - O , (4) 

1 k 1 



with Gjj(k) = —T l 6&/6"Lji(k). Equation (ffj) corresponds to 

©[£] = - \ y w(k - k'm [t(k)p 3 t(k')P3] , (5) 

kk> 

where W(k - k') = Yjx V~ 1 (x)e ~ !( * is the matrix inverse 
of V(k - k') in the (k, k') space (W > represents an attrac- 
tion). The stationary condition for £2[£] again gives the Dyson 
(Eliashberg) equation. Since the anomalous component of the 
self-energy matrix is the superconducting order parameter, (O 
is regarded as the Landau free-energy functional for supercon- 
ductivity. It is straightforward to extend the present formula- 
tion for non-scalar general interaction by modifying only the 
part of ©[£] in the free-energy functional, Q[S], 

Once the stationary solution of is obtained, the physical 
equilibrium free energy is given by 



n. = -\ Tr ln H^W) + \c{k)t(k) 



- Oo, (6) 



in which the Dyson equation, G~ l (k) = G Q l (k)-'Z(k), is satis- 
fied. It is useful to note that the derivative of Q. s with respect 
to a parameter x can be obtained by the explicit derivative of 
Q. evaluated at the stationary, i.e., dQ. s /dx = dtl/dx^ . 



B. Decomposition of singlet and triplet components 

Let us introduce the components of the 4x4 matrix as fol- 
lows, 

r-'m ton - ( z(k) - A(k) \ 

z(k) = zo(k)cr + z(k) ■ a- = -z + (-k) T , 

zo(k) = ico n -%k- Zo(*), z(k) = h - £(*), 
A(k) = do(k)T Q + d(k) ■ t = -A(-k) T , 
A + (k) = d^(k)T f Q + d + (k) ■ r f = -A + (-k) T , (7) 

6(k) = -(T T ^(k)) = (_f^ k) 

G(k) = go(k)cr + g(k) ■ <x = -G + (-k) T , 
F(k) = Mk)r + /(*) ■ T = -F(-k) T , 



F + (k) = f+(k)Tl + f + (k) ■ r t = -F + (-k) 7 



(8) 



where the superscripts f and T represent the hermite conju- 
gate and the transpose of a 2 x 2 matrix, and and h are 
the one-particle energy measured from the chemical poten- 
tial and the external magnetic field, respectively. The normal 
part is decomposed into the charge and the spin components 
by <t and cr, which are the 2x2 unit matrix and the vector 
of the Pauli matrix acting on the spin space. The relations, 
z (k) = -z + (-k) T and G(k) = -G + (-k) T are obtained by defi- 
nition of the Green's functions. The anomalous self energy is 
also decomposed into the singlet and the triplet components 
by to = icTQ(T2 and r = !<to"2 . The anti-commutation relation 
of fermions requires the even property, do(k) = do(-k) and 
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d^ik) — d^i—k) in the singlet channel, and the odd property, 
d(k) = -d(-k) and d + (k) = -d + (-k) in the triplet channel. In 
this paper, we assume the presence of the inversion symme- 
try, and z + (k) = -z(k)\ i.e., z + a (k) = -z a (k) m (a = 0, 1,2,3) 
holds. Note that the d^(k) [/+(£)] are independent of d a (k) 
lfa(k)], and the relation between them are determined by real 
and minimum condition of the free-energy functional as fol- 
lows. 

With this preliminary, we expand the free-energy functional 
(HJi with respect to do(k), d(k), d+(k), and d + (k). After some 
manipulation (see, Appendix in detail), we obtain the lowest- 
order expression of the Landau expansion as 

An 2 [d a , d+] = w ( k ~ k ') [do(k)d^(k') + d(k) ■ d + (k')] 

kk' 

+ tY j -^[(zozS - z ■ z + )(ckdl + d-d + ) 



+ (z x d) ■ (z + x d + ) + (z x d + ) ■ (z + x d) 



no ■ mo + n\ ■ m\ + n 2 ■ m 2 



(9) 



where the free-energy functional is measured from that of the 
normal state, and w - \z\, z a and z„ are evaluated in the normal 
state. The symmetric and the antisymmetric vectors have been 
introduced as 



no(k) = zo(k)z + (k) + z + {k)z{k) = n (-k), 
ni(k) = i[z(k) x z + (k)] = -ni(-k), 
n 2 (k) = z (k)z + (k) - z + o(k)z(k) = -n 2 (-k). 



(10) 



Due to the relation, z^(k) = -z a (k)*, no(k) and n\{k) are real, 
while n 2 (k) is pure imaginary. The corresponding symmetric 
and antisymmetric vectors composed of d a (k) and d^(k) are 
given by 



mo(k) = i[d(k) x d + {k)] = mo(-k), 

mi(k) = dl{k)d(k) - do(k)d + {k) = -mi(-k), 

m 2 (k) = dl{k)d(k) + do(k)d + {k) = -m 2 (-k). 

The explicit expressions of n,(fc) are given by 



1 



n (k) =%h + T.' Q h - ^S' + w„S" - (E S*)', 



m(k) = 2(h x S") - i(S x S*), 
1 



(11) 



(12) 
(13) 



-n 2 (k) = -oj n h + 2 ( > + ^S" + co„S' - (EoS*)", (14) 

where the prime and the double prime mean the real and the 
imaginary part, respectively. Since the normal self energy 
does not change the symmetry of the system in the normal 
phase (we do not consider a spontaneous symmetry breaking 
in the normal self energy), we neglect the normal self energy 
in considering symmetry properties of m(k). Then, no ~ %h, 
n\ ~ 0, n 2 ~ ico„h, and if one of them is finite, a mixing 
between different components of d a arises. 

Let us first consider the case h — 0, i.e., n,(Af) = 0. In 
this case, either the pure singlet or the triplet pairing is real- 
ized depending on the structure of the attraction. Since the 



equilibrium free energy must be a real quantity, a stationary 
solution must have the relation, d^{k) = (f> a d a (k)* with (p a tak- 
ing either +1 or - 1 . Note that (f> a should be independent of k, 
since d a (k) and d^{k) are mutually equal counterparts. For a 
pure singlet or a triplet state, it is the sufficient condition for 
the real free energy. Then, the sign <p a = + 1 should be cho- 
sen by the minimum condition of the free-energy functional 
through a phase transition to the superconducting state. Con- 
sequently, da(k) = d a (k)* always holds for the pure singlet 
or the triplet pairing, irrespective of its symmetry in the fre- 
quency domain, as was discussed previously 19 . Hereafter, the 
free-energy functional, that is defined over the hyperplane un- 
der the constraint between d a (k) and d^(k), is referred as the 
constrained free-energy functional. Note that the constrained 
free-energy functional is always real by definition. 

Next, we elucidate the effect of the external magnetic field, 
h. In this case, Re(n-o) and Im(n2) become finite (no = in 
the presence of the particle-hole symmetry). Thus, the terms 
coupled with no and n 2 in the free-energy functional induce 
mo and m 2 . Due to the real condition of the free energy, 
mo and m 2 must be real and pure imaginary, respectively. 
In order to satisfy these conditions, the singlet and the triplet 
components have the opposite sign relation, i.e., 

d^{k)=<pd {kT, d + {k) = -<f>d(k)\ (0 = +1, or - 1). 

(15) 

It should be noted that n 2 does not mix the spatial parity, but 
mix the even-frequency and the odd-frequency components. 
Namely, either the singlet or the triplet has the odd-frequency 
dependence. By inserting (TT3T > into ( fTTT i. it is easily shown 
that the relative phase factor between the singlet and the triplet 
components is pure imaginary. 

Using ( fT5l ), the explicit form of the equilibrium free energy 
is eventually given by 

£2 S = - T - 2 InX(fc) - | Ylgo(k)Z (k) + g(k) ■ S(fe) 

k k 



- 4>{Mk)d (kT - f(k) ■ d(ky 



-n , (16) 



where g a (k), I< a (k), f a (k) and d a (k) satisfy the Eliashberg 
equations, ( IA12I ), ( IA13b with ( IA24I ) and ( IA25I ). The explicit 
expression of X(k) = \-G~ l (k)\ is given by dAlOt . The anoma- 
lous Green's functions satisfy the relations similar to (fT3T > as 



fo + (k) = Mm*, rm = -*/(*)•. 



(17) 



The sign of <p is determined by the minimum condition 
of the constrained free-energy functional, or equivalently the 
solvability of the Eliashberg equations in the superconducting 
state. As will be shown explicitly in the next section, a sta- 
ble solution of the Eliashberg equations can be found, only 
when we use the correct choice of the sign, <j>. Note also that 
this stable solution corresponds to the saddle point of the con- 
strained free-energy functional. An arbitrary initial condition 
always converges to this saddle point in the iterative proce- 
dure, as long as the correct choice of (f> is used. In a coex- 
istence case and one component dominates over other com- 
ponents, the majority component would determine the sign of 
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<f>, i.e., <p = +1 for the singlet majority, while <f> — -1 for 
the triplet majority. In the competing case, however, the com- 
parison of the equilibrium free energies both for <p = ±1 is 
necessary to determine the correct sign of <p. 

It is important to note that the contributions to the equilib- 
rium free energy from the singlet and the triplet components 
have opposite signs as in the last term of ( fToT l, as a conse- 
quence of (fT5l l. In other words, the mixing of the singlet 
and the triplet pairings tends to increase the net free energy 
as compared with the case of neglecting the odd-frequency 
component. As a result, it alters behaviors of various physical 
quantities such as the reduction of T c . The explicit examples 
are given in the next section. 



III. THE v-WAVE SINGLET WITH ODD-FREQUENCY 
TRIPLET STATE UNDER MAGNETIC FIELDS 

In order to elucidate essential features of the even- and the 
odd-frequency mixing, we consider the s-wave singlet and 
triplet pairings under magnetic field h along z axis as the sim- 
plest example. The orbital effect is not taken into account for 
simplicity, so that the discussions can be applied in low fields. 

A. The formulation based on the free-energy functional 

Since «2 C ~ iu>„h + in this case, the singlet (do) and the 
S z = triplet (dj) components are mixed 26 , and the relative 
phase becomes pure imaginary as noted in the previous sec- 
tion. We assume the particle-hole symmetry after perform- 
ing the momentum integration, the no term in the free-energy 
functional vanishes. By these reasons, we denote do(icj„) and 
di(ia>n) as d s „ and id t „, then both d s „ and d ln can be chosen as 
real. The anti-commutation relation requires that d s „ and d tn 
have the even-frequency and the odd-frequency dependences, 
respectively. Therefore, without considering the retardation 
effect (neglecting aj„ dependence), d m should vanish as in past 
studies. 

The normal part of the self energies T.o(icj„) and E3(;w„) are 
finite as well. It is shown that Eo(ia>„) [I,3(iu>„)] is pure imag- 
inary odd [real even] function, so that we denote Eo(ia) n ) = 
ito„{\ - Z„) and E3(/«„) = £3,,. Then, both Z„ and l,^ n are 
the real even functions. Z„ is so-called the mass enhancement 
factor. 

For notational simplicity, we introduce A„ = d m , + id^, = 
A* _j = <pA^ corresponding to the (T,i) component, and 
E„ = iu>„(l - Z„) + £3,, = S* _j corresponding to the (T,T) 
component. We adopt the local electron-phonon-type attrac- 
tion with retardation, 



PfV„ ui 
A 



— (W )m,« — 



Wg + (u>„ 



0J„) Z 



(18) 



where u F is the characteristic range of the interaction, p F is 
the density of states per spin at the Fermi energy, and A > is 
the dimensionless coupling constant for the attraction. 

The Einstein-phonon attraction ( fT8l is the same order of 
magnitude A both in the even-frequency (the even part with 



respect to to,, — > -w„) and the odd-frequency (the odd part 
with respect to to„ — > -w„) channels (hence the same order of 
T c ), unless the normal self-energy effect (mass enhancement) 
is not taken into account 12 . Owing to the odd property of the 
attraction in the odd-frequency channel, which is vanishing 
toward T — > 0, the odd-frequency pairing is suppressed at 
low temperatures, yielding the reentrant behavior in the case 
of the pure odd-frequency pairing. In addition to this, the 
mass enhancement significantly suppresses the odd-frequency 
pairing, so that it is naturally expected that the induced odd- 
frequency component discussed below is also small as com- 
pared with the majority even-frequency component^. 

Performing the k integration, we obtain the free-energy 
functional as 



= T E Wm '" ( ^ AmA " _ * m * n) ~ 2 X X ln x(k) ~ Q ° 

mn n h 

= -jX W ">>» (0 A m A « - £ m £„) - 2p F nT ^(A, - \cj n \), 

mn n 

(19) 



where X(k) = |£ 2 + D 2 | 2 with 



D„ = yj(oj„ + i'E„ - ih) 2 + 0A 2 = D*_ n _ v 



(20) 



The stationary condition, 6Q./61,,, = 6Q./6A,, = leads to 
the Eliashberg equations, 



AnTY 



A m = AjtT ^ v m , n F n , 

n 

where the local Green's functions are defined as 
~1 d . „,.. _lo„ + iL„-ih 



(21) 



-1 v o io n + il„ — in 

2np F ^ 3S„ D„ 



2np F dA n D„ 



(22) 
(23) 



These Eliashberg equations are examined numerically in the 
previous work 1 -, Using the solution of the Eliashberg equa- 
tions, the equilibrium free-energy difference is expressed as 

AQ S = 2p F nT V Re [(cf>F n A n - G„S„) + w„(l - Z N „) 



n=0 



-2(D„ - Zn„w„)] , 



where Z Nm = 1 + (AnT/\u> m \) Z7=o( v m,n - v m ,-„_i) is the mass- 
enhancement factor in the normal state. By appropriate nu- 
merical differentiation of AQ S , we obtain the thermodynamic 
quantities such as the specific heat. The superfluid density 
may be obtained by the standard derivation a s 19 ' 23 



« s = 2*rf>(^L 



n=0 

OO 

=^E 

n=0 



(24) 
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which can be applied in weak magnetic fields, since we have 
neglected orbital effects. There exist the paramagnetic and 
the interference contributions to n s other than the ordinary dia- 
magnetic one, since the signs of d% and d\ terms are opposite. 

In what follows, we adopt relatively large A in order to em- 
phasize characteristic features of the even-odd-frequency mix- 
ing. However, qualitative features of the results are not altered 
for weaker A and h. We use the unit of energy aswe = 1. 



B. Single-gap approximation in the strong-coupling limit 

Let us first consider the strong-coupling limit, A — > oo in or- 
der to grasp an outline of the solutions. In this limit, n — and 
-1 components, i.e., Ao = A*j eA s + ;A t dominate over other 
components 27,28 . As will be shown below, this simple approx- 
imation provides important information on the structure of the 
free-energy functional and qualitatively similar behaviors to 
those obtained by solving the full Eliashberg equations 1 -. The 
approximate free-energy difference is given by 



0.02 



2p F /A 2 A 2 \ 
AQ[A s ,A t ] = — - — \-4p F co Re[D -Zo<o ]. (25) 

Here, the self energy is evaluated in the normal state for sim- 
plicity as Z = 1 + A. D - ^(Zqojq - ih) 2 + 0(A S + /A t ) 2 . 
The normalized attractions in the singlet and the triplet chan- 
nels are introduced as 



U s (t) = V ,0 + V -1 = 



«t(0 = vo,o - v o,-i 



2 + f 2 
1 + t 2 ' 
t 2 

l+t 2 ' 



(26) 
(27) 



with t = 2jtT/coe- The corresponding gap equations and the 
equilibrium free-energy difference are given by 



A s = Au s u>{) Re 
A t = Auitoo Im 



A s + /A t 



A s + /A, 



AQ S = 2p F w Re 



Do 

</>(A s + ;A,) 2 - w (l - Zo)(Z w - ih) 



(28) 



Do 

+ a>o(l - Zq) - 2(D - Zqcjo) 



(29) 

Expanding Q with respect to A s and A t , we obtain the 
lowest-order Landau free energy as 



AQ[A s ,A t ] = 2p F 0(A s , A t )7C^J+0(A 4 ), 



r K{T,h) 



l/Au s - A 
B 



with 



Z w ( 2 



(Zou ) 2 + h 2 - 



B 



B 

(1/Aut-A) 



u>oh 
(Zo(Oo) 2 + h 2 ' 



(30) 
(31) 

(32) 
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FIG. 1 : (Color online) The T dependence of eigenvalues of 'KiT, h). 
One eigenvalue is always negative, while the other changes its sign 
at T c (h) except for h/a>E = 3. 



Note that the minus sign appears in the (2,2) component of 
PTI ) due to the relation, $15[ . The superconducting instabil- 
ity at T c {h) is determined by the condition |7C(r c , h)\ = 0, or 
equivalently, by the condition such that one of the eigenvalues 
becomes zero. 

FigureQ]shows the T dependence of eigenvalues of 'KiJ, h) 
at Ii/u>e — 0, 1,2, and 3 for A = 10. One eigenvalue coming 
mainly from the triplet channel is always negative as a con- 
sequence of ( fTBI ). while the other changes its sign at T c (h). 
This behavior indicates that the normal state is not character- 
ized by the minimum as usual, but by the saddle point of the 
constrained free-energy functional above T c . Then, the sad- 
dle point at the origin turns into the unstable maximum (stable 
minimum) below T c in the case of = +1 ((f) — -1). There- 
fore, the choice of <\> — + 1 describes a realistic phase transition 
to the superconducting state. In other words, (f> — +1 is cho- 
sen since the singlet pairing is the majority component in the 
present case (the pure triplet pairing is suppressed completely 
by the mass-enhancement effect). In the absence of the mag- 
netic field where the singlet and the triplet channels are de- 
coupled, the sign (f> can be chosen independently in the singlet 
((f) = +1) and the triplet ((f) = -1) channels. However, the 
finite coupling in the presence of the magnetic field requires 
the opposite sign of </> between two components, yielding the 
energy loss in the triplet channel, even though the interac- 
tion is originally attractive in the triplet channel. Although 
the relative-sign requirement ([T5t leads to higher free energy 
than the case of neglecting the odd-frequency triplet compo- 
nent as in the past studies, A t = line in A s -A t plane is no 
longer a stationary solution of the free energy in the presence 
of magnetic fields, and the saddle point of the constrained 
free-energy functional gives a stable solution. The upturn of 
the eigenvalue in lower temperatures suggests a reentrant be- 
havior of the superconducting phase. For higher fields (e.g. 
h = 3), no sign changes occur and the system remains the 
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normal state. 
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FIG. 2: (Color online) The T dependence of two real eigenvalues ( 
of the kernel in the linearized gap equation. The dashed line with 
the label h c represents £ for h c (T) where two eigenvalues coincide. 
((T c ) = 1 gives the transition temperature. 
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which are shown in Fig. [2] for /i/we = 0, 1,2 and 3 at A = 
10. The dashed line with the label h c represents £ for h c (T) 
where two eigenvalues coincide. In the presence of magnetic 
fields, the kernel of the linearized gap equation becomes non- 
hermite, leading to the complex eigenvalues when 



h> h c = 



*V2 



(35) 



+ t L 



The eigenvalues shown in Fig. [2] provide qualitatively simi- 
lar behaviors to those obtained by solving the full linearized 
Eliashberg equations, (l2"TU^. 

The T-h phase diagram determined by \K(Tc, h)\ — or 
£(T ct h) = 1 for A = 5, 10, and 15 is shown in Fig. [3] As 
was expected, the superconducting phase boundary shows the 
reentrant behavior in high fields. T c (li) is considerably sup- 
pressed as compared with the case of neglecting the odd- 
frequency triplet component. 

Figure [4] shows the typical landscapes of the constrained 
free-energy functional AO[A s , A t ]/2pp for A = 10. The upper 
panels (a)-(c) show the T dependence at h/a>E = 0.5, and the 
lower panels (d), (e) show the h dependence at T/u>e = 0.1. 
FigureUff) is the case of <p = — 1 at T/we = 0.1 and /i/o»e = 
1 .5. The closed circles represent the saddle points correspond- 
ing to the solutions of the gap equation (see, Fig. [5]). The open 
circle and squares indicate other saddle point and maxima, 
whose free energies are higher than the saddle points repre- 
sented by the closed circles. The closed square represents the 
minimum of Afi in the case of (/> = — 1. By the compari- 
son between Figs. [4] and |5j we can see that both the normal 
and the superconducting states are characterized by the sad- 
dle points of the constrained free-energy functional with the 
choice of <p = +1. Note that an arbitrary initial condition 
always converges to the solution of the Eliashberg equations 
that corresponds to the saddle point of the constrained free- 
energy functional in the case of <p - + 1 . On the other hand, 
in the case of (p = — 1, the energy of the saddle point is always 
higher than that of the origin in the constrained free-energy 
functional, which corresponds to the normal state. This is the 
consequence of the wrong choice of (p. 



C. Comparison with or without the triplet component 



FIG. 3: (Color online) The T-h phase diagram for several coupling 
constants. The phase boundary shows reentrant behavior in high 
fields. 

The superconducting instability can also be observed by the 
linearized gap equation, 



/A s \ _ u s A -u s B\ /A s 

Ma, \u t B Ut A U t r 



(33) 



where the eigenvalue £ = 1 at T c signals the superconducting 
instability. The eigenvalues are given by 



£ = AA 



1 f 2 (2 + t 2 )h 2 

1 + z k Jl 



1 + f 2 



(Z w ) 2 



(34) 



Let us now discuss some physical quantities by solving the 
full Eliashberg equations, ( 1211 1. and compare with those ob- 
tained by neglecting the odd-frequency component (Im A„ = 
d tn is forced to be zero during the computation), which corre- 
sponds to the limit of vanishing attraction in the triplet chan- 
nels. Figure|6]shows the comparison of the specific heat C(T) 
and the entropy difference AS (T) from the normal state. Here, 
7n = 2pF^ 2 (l + A)/3 is the renormalized Sommerfeld coeffi- 
cient in the normal state. Although T c (h) is considerably dif- 
ferent from each other, the discontinuities at T c show almost 
no difference {AC/jnT c ~ 1) and the slopes of C(T) below T c 
differ slightly from each other. Namely, the entropy release in 
the mixed-frequency solution is more gradual than that in the 
pure singlet solution. Note that the correct mixed-frequency 
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(a) T=0.1, h=0.5 



(b) T=0.3, h=0.5 



(c) T=0.5, h=0.5 




FIG. 4: (Color online) The landscape of the constrained free-energy functional in unit of 2pp measured from that of the normal state, (a)-(c) 
the T dependence at h/cj E = 0.5, (d), (e) the h dependence at T/o>e = 0.1, (f) the case of <j> = -1 at T/u> E = 0.1, h/a> E = 1.5. The closed circles 
represent the saddle points corresponding to the solutions of the gap equation (see, Fig. [5}. The open circle and squares indicate other saddle 
point and maxima, whose free energies are higher than the saddle points indicated by the closed circles. The closed square represents the 
minimum of AC1 in the case of <p = - 1 . On the white line in the right-lower region in (f), Re(Z)o) = at which the derivative is discontinuous. 
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The comparison of the superfiuid density n s (T) is shown 
in Fig. [7] The evolutions of n s (T) near T c are linear in both 
cases, but the slope of the mixed-frequency solution is slightly 
smaller than that of the pure singlet solution. In low enough 
temperatures, however, the superfiuid density of the mixed- 
frequency solution is larger than that of the pure singlet so- 
lution. This behavior can be interpreted by a paramagnetic 
contribution from the induced odd-frequency component in 
the vicinity of T c , which suppresses the Meissner screening 
slightly. The majority singlet component gives the ordinary 
diamagnetic Meissner current, which eventually dominates 
over the paramagnetic contribution from the induced odd- 
frequency component at low temperatures. The existence of 
the paramagnetic Meissner effect is extensively discussed in 
the spatially nonuniform systems at the surface/interface^^. 



FIG. 5: (Color online) The T dependence of the gap functions and 
the equilibrium free energy in unit of 2p F at h/a) E = 0.5. 



solution has lower entropy than that of the pure singlet solu- 
tion. 



IV. SUMMARY 

We have investigated general structure of the free-energy 
functional for the coexistence state of the singlet and triplet 
pairings with the even-odd-frequency mixing. The real con- 
dition of the free-energy functional for the coexistence state 
requires the relation between the gap functions and their 
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free-energy functional. The minority component with the ex- 
traordinary relation, d^(k) = -d a (k)* suppresses the super- 
conducting state yielding lower T c as compared with those 
obtained by neglecting the minority component. It also re- 
flects in property of various physical quantities such as the 
specific heat and the superfluid density. 

Using the simplest example of the ordinary s-wave singlet 
state mixed with the triplet state with the odd-frequency de- 
pendence under magnetic fields, we have demonstrated the 
saddle-point structure of the constrained free-energy func- 
tional, and the induced odd-frequency triplet component (the 
minority component with the opposite sign) alters the T de- 
pendences of the specific heat and the superfluid density 
as compared with those obtained by neglecting the odd- 
frequency triplet component. It would be interesting to de- 
tect such effect under the broken time-reversal symmetry in 
strong-coupling superconductors. 



FIG. 6: (Color online) The comparison of the /"-dependences of the 
specific heat with or without the odd-frequency component. The en- 
tropy difference is shown in the inset. 



Acknowledgments 



0.1 



0.05 



I 


I I I 

— even-co + odd-(o 




— even-co only 


- \ = 10, h=1 







1 1 1 i 1 1 1 i 1 1 i 1 1 i 1 1 1 i i 1 1 i 1 1 1 1 1 MM 1 1 i 1 1 i 1 1 1 



0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 
T/T c 

FIG. 7: (Color online) The comparison of the Independence of the 
superfluid density with or without the odd-frequency component. 



We acknowledge valuable discussions with Y. Fuseya, K. 
Miyake, and Y. Tanaka. This work is supported by a Grant-in- 
Aid for Scientific Research C (No. 23540414) from the Japan 
Society for the Promotion of Science. One of the authors 
(H.K.) is supported by a Grant-in-Aid for Scientific Research 
on Innovative Areas "Heavy Electrons" (No. 20102008) of 
The Ministry of Education, Culture, Sports, Science, and 
Technology (MEXT), Japan. 



Appendix A: Derivation of the free-energy functional, and the 
Green's functions 

1. The free-energy functional 

Here, we show the derivation of the free-energy functional 
in terms of the spin-decomposed components. Let us start 
from the expression (0). For ©[£], it is straightforward to ob- 
tain, 



particle-hole converted counterparts as, 

dJ(Jfc) = (f>do(k)*, d + (k) = -4>d(kY, 

where <p takes either + 1 or - 1 depending on which component 
is majority, i.e., <p = +1 for the singlet majority and (f> = -1 
for the triplet majority, in order to lower the resultant equilib- 
rium free energy via the interference effect. In the pure singlet 
or the triplet state, this requirement reduces to the ordinary 
relations^, d+(Jfc) = d (k)* or d + (k) = d(k)*. 

The opposite signs in the singlet and the triplet channels 
lead to the anomalous structure of the free-energy functional, 
which is constrained by the above relations, and the equilib- 
rium normal and superconducting states are characterized not 
by the minimum, but by the saddle point of the constrained 



kk' 

4Zj W(k K)lT [\A + (k) Z + (k))\-A + (k') Z + (k')J 

= _I V W(k _ k > )tx [X(fc)E(*0 - A(k)A + (k')] 

1 kk' 

= Y J W(k- k') [d (k)d^(k') + d(k) ■ d + (k') 

kk' 

-EoTO(fc') - £(*) • £(£')] ■ (Al) 

In the first term of ©, we use the identity, Trln(- • ■ ) = 
lndet(- ■ ■ ). For a general 4x4 matrix that consists of 2 x 2 
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matrices, A, ■ ■ ■ ,D, can be decomposed as 

I A B\_ll BD l \lA{l-A- x BD l C) 0\[ 1 
C D J~ 1 D fl 'C 1 



(A2) 

Using this identity, and omitting k for notational simplicity, 
we obtain 







(a* -z\ 





i - (zr 1 A(z + r l A + 



where 



w = \z\ = zl - z\ w + = |z + | = zf - z +2 . 



By the straightforward calculation, we have 

(z) _1 A = — [or T + a-r], 
w 

^" lA+ss ^K T S +a+ - Tt ]- 

where we have defined 

ao = zodo - z ■ d, 

a = zod - zdo - i(z x d), 

a o = z o d o + z+ • d+ , 

a + = + z + d + Q + x d + ). 

Using these expressions, we obtain 

Dqctq -{D + M) a 



(A3) 

(A4) 

(A5) 
(A6) 



p = i -(zr 1 A(z + r'A + 

where 



Dq = ww + - (ao^o + a ' a+ )' 

D = q'oq; + + aaj, A-f = /(a x a + ). 



(A7) 



(A8) 



(A9) 



Putting these expressions into dA3t and using D ■ M = 0, we 
have 



X(k) = 



Dl-D 2 - M 2 



(A10) 



The final expression of the free-energy functional is given by 
Q[£] = V#- fc') [j (^K(^') + d(k) ■ d + (k') 

kk' 

-lo(fcW) - m) ■ - | J] In , (Al 1) 

where Xo(k) is the non-interacting value of X(k). 

The stationary condition of £2[£] gives the Eliashberg equa- 
tions, 

Z a (k) = ?rJ]V(k-k')g a (k'), (a = 0,1,2,3), (A12) 

1 k' 

d a (k) = ?-Vv(k-k')f a (k'), 

1 k> 

<£(*) = \ Yj V(k ~ ^ti&l (A13) 



where we have defined the Green's functions as 
S]nX(k) 



ga(k) = - 

fjk) = 



SZaik) 

6\nX(k) 



Sd+{k) 



dlnX(k) 
6d a (k) 



(A14) 
(A15) 



Note that 'E(k) and X + (£) in X(k) should be treated as inde- 
pendent quantities in the functional derivative. Their explicit 
expressions are given later. 

Using the solution of the Eliashberg equations, we obtain 
the equilibrium value as 

I4z [/o(*K(*) + /® • d+ w 
z 1 k 

- g0 (k)Z (k)-g(k)--£(k)], (A16) 

and the explicit equilibrium free energy is given by (ITol with 
use of the relation, (1151 . 

As we discussed in detail in §11. B, the real condition of the 
free-energy functional requires the relation, (Q3}, 

do(k) = <f>do(kT, d + (k) = -<pd(k)*, (<f> = +1, or - 1). 

With these constraints, we have 

a o — -<P a {)> a+ - <P a *- (A17) 
From (IA13b and (fl3T l. it can be shown the relation, (flTt . 
f Q + (k) = <PMkf, f + (k) = -</>f(k)*. 

2. The lowest-order Landau expansion 

In order to obtain the lowest-order expression of the Landau 
expansion, we expand the last term of ( IA1 U with respect to 
a a and a* (a — 0, 1 , 2, 3), then we have 

--In 



2 \X (k)j \ww + ) 



= -Tin I I - 
T 



a a 



+ a a 



(aoffo + a ' a+ ) > (A18) 



where cn = -(T/2) ln(vvvv + /^o)- Using the explicit ex- 
pression, 

OoOq +aa + - (zoZq - z ■ z + ){d {) dl + d ■ d + ) 

+ (zxd)- {z + x d + ) + {zx d + ) ■ (z + x d) 

- no • mo + n\ ■ m\ + n 2 ■ rri2, 

= (zoZq - z ■ z + )d du + (zoZq + z ■ z + )d ■ d + 

j 1,2,3 

- 2 2 Ni J Mi J ~ n ° ' m ° + n i ' m i + n 2 ■ m 2 , 

(A19) 

where n,- and m, are defined by ( fTOb and (fTTT i, we obtain the 
lowest-order Landau expansion as ©. The symmetric matri- 
ces are defined by Ny = z,z| + z\zj and My = + d^dy, 
respectively. 
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3. The Green's functions 

The stationary condition of the free-energy functional leads 
to the Dyson (Eliashberg) equation, 



G{k) -F(k)\_( z(k) -A(k)Y 
-F + (k) G + (k))-\-A + (k) z + (k)j 



(A20) 



Thus, by inverting the matrix explicitly, we obtain the expres- 
sions of the Green's function. For this purpose, we again use 
the identity, 1A2) . then we express the inverse of a matrix as 



1 0\ 1 (A-BD l C 0\ '/l BD 1 ^ ' 



D l C 1 







D 1 



1 0\((A- BD l C) x \(\ -BD 

1 



-zr'c ij\ o d 'Ho 

(A - BD l C) 1 {C-DB x A) l \ 
(B-AC l D) 1 (D-CA- l B)- l > 



Using this formula, we have 



(A21) 



(A22) 



G(k) = [l - (zr 1 Mz + r 1 A + l (z)- 1 = (zP)-\ 

F(k) = -[l- (zr'A^T'A^ ' (z)" 1 ^*) -1 = -GA(z + r'. 

(A23) 

Using dA8b and after some manipulations, the explicit expres- 
sions are given by 

go(k) = ~ [w + zo - (dod^ + d ■ <i + )zo + (mo - m 2 ) ■ 2 + ] , 



-w + 2; + (dodg - d ■ d + )z + + (mo + m2)zo 



+ (z + • d)d + + (z + ■ d + )d + i(mi x 2 + ) 



(A24) 



-{d 2 - d 2 )dl + (zoZq - z ■ z + )d + (ni + n 2 ) ■ d 
(dl - d 2 )d + + (zoZq + z ■ z + )d - (ni - n 2 )do 
(z + ■ d)z - (z ■ d)z + - /(n x d)\ (A25) 



Appendix B: Effect of the Rashba spin-orbit coupling 

Here, we consider briefly the effect of the Rashba antisym- 
metric spin-orbit coupling 16 , jk = -~f-k = An(k x z) — 



Asiky, —k x , 0). In this case, we just replace the spin part with 

z(k) = - 7fe - E(*). (Bl) 

In the presence of the Rashba spin-orbit coupling, z + (k) = 
-z(k)* does not hold in general, but ww + = \w\ 2 still holds. 

The symmetry properties of the symmetric matrix 2Vy and 
the symmetry-lowering vectors n, in ( |A19t are extracted as 



N(k) ~ Al 



&y k x ky 

0, 
n (k) ~ iaj„j k = ;'w„/l R (^, ~k x , 0), 
n x (k) ~ 0, 

n 2 (k) ~ %klk = tk^Riky, -k x , 0). 



(B2) 

(B3) 
(B4) 
(B5) 



In order to gain the interference energies, and the free energy 
is to be real, we should have the relation, 



<% 12 (k) = <t>d oxz (ky, dim = -<t>d 3 (ky. 



Then, the conjugate quantities become 



(B6) 



M(k) = 2<p 



Rs(did~) -ilm(d\dX) 

" 2 



Reid^l) \d 2 \ -ilm(d 2 dl) 
-i\m{d\d^) -ilm(d 2 d^) -\d^\ 2 

mo(k) - -2(p[iR&{d 2 dl),-iRe(d2d\),\m(did*2)], 

m\{k) = -2(f>[ilm(dodl),ilm(d()dl),-Re(dQd^)], 

m 2 (k) = 20[Re(c/o^iX Re^oa^X ~ l \vnid a d*^)\. 



(B7) 

(B8) 
(B9) 
(B10) 



The relative phases among d a should be taken as real, and 
lm(d a dp vanishes. Note that the frequency mixing occurs be- 
tween (do,di,d 2 ) and di components via the no term, while 
the spatial parity mixing occurs between (c/0,^3) and (d\,d 2 ) 
components via the Nij, no and n 2 terms. Consequently, the 
spin mixing occurs between do and d components. For in- 
stance, it is possible to emerge the coexistence such as 

do : an s-wave singlet with the even-w dependence, 

d\ : a p-wave triplet with the even-w dependence, 

d 2 : a /?-wave triplet with the even-w dependence, 

dj : an s-wave triplet with the odd-w dependence, (Bl 1) 

purely by the symmetry consideration. 
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